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Abstract 

Following Waxman and Bahcall [Q] we calculate the event rate, energy and zenith 
angle dependence of neutrinos produced in the fireball model of gamma ray bursts 
(GRB). We emphasize the primary importance of i) burst-to-burst fluctuations and 
ii) absorption of the neutrinos in the Earth. From the astronomical point of view, 
we draw attention to the sensitivity of neutrino measurements to the boost Lorentz 
factor of the fireball T, which is central to the fireball model, and only indirectly 
determined by follow-up observations. Fluctuations result in single bursts emitting 
multiple neutrinos, making it possible to determine the flavor composition of a beam 
observed after a baseline of thousands of Megaparsecs. 

PACS number(s): 96.40.Tv, 98.70.Rz, 98.70.Sa 
Keywords: Gamma Ray Bursts, Neutrino detection 

1 Introduction 

The origin of GRB is one of the most fascinating outstanding problems in astronomy. Their 
observed energy injection in the Universe is sufficiently large to possibly resolve another 
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long-standing puzzle: the origin of the highest energy cosmic rays 0, |3| . Mounting evidence 
suggests that GRB emission is produced by a relativistically expanding fireball, energized 
by a process involving neutron stars or black holes ||. In the early stages the fireball, its 
radiation trapped by the very large optical depth, cannot emit photons efficiently The 
fireball's kinetic energy is therefore dissipated until it becomes optically thin — a scenario 
that can explain the observed energy and time scales of GRBs, provided the bulk Lorentz 
factor of the expanding flow, T, is > 10 2 — 10 3 . 

Protons accelerated in shocks in the expanding fireball interact with photons to produce 
charged pions, the parents of high-energy neutrinos [p], ||- Assuming that particles acceler- 
ated in the GRB sources produce the observed cosmic rays above the "ankle" of the energy 
spectrum near 3 x 10 18 eV, one derives that the average single burst produces only ~ 10~ 2 
neutrino events in a high energy neutrino telescope with 1 km 2 effective area. Although 
the expected rate is therefore low, the neutrino signal should be relatively easy to observe 
provided the detector is large enough: GRB neutrinos will have a hard spectrum extend- 
ing well beyond the background from atmospheric neutrinos and, even more important, the 
high-energy GRB neutrino events should coincide with observed GRB photon events within 
a narrow time window. 

In this paper we calculate the experimental signatures of GRB in a kilometer-scale neu- 
trino detector such as the proposed IceCube 0. We emphasize the importance of taking 
into account burst-to-burst fluctuations [0] as well as absorption of the neutrino signal in 
the Earth for both event rates and experimental signatures. Both effects produce additional 
and striking signatures with discriminating sensitivity to the value of the bulk Lorentz factor 
whose value is only indirectly inferred from other astronomical observations |fj, |S| . 

The observation of GRB neutrinos over a cosmological baseline has scientific potential 
beyond testing the "best-buy" fireball model: the observations can test with unmatched pre- 
cision special relativity and the equivalence principle, and study oscillating neutrino flavors 
over the ultimate baseline of z ~ 1 [p]]. 
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2 Calculation of GRB Neutrino Rates and Signatures 



In calculating the event rates and experimental signatures of GRB neutrinos in a high en- 
ergy neutrino telescope we follow the model of Waxman and Bahcall |l] as implemented by 
Halzen and Hooper [Tj. We have normalized the neutrino flux to the energy rate injected 
in the Universe needed to explain the observed cosmic ray (CR) spectrum above 10 19 eV, 
Ecu = 4 x 10 44 ergs Mpc -3 yr~ 4 . This energy rate was calculated in reference ||, assuming 
a cosmological distribution of sources and taking into account CR propagation in the Cosmic 
Microwave Background Radiation (CMBR). The value quoted above corresponds to the "low 
redshift" (z < 1) energy generation rate of CRs. Note that, because of the absorption of 
highest energy cosmic rays by CMBR photons, one could further increase Eqr without di- 
rectly affecting their observed flux. Waxman and Bahcall [10] have calculated an upper limit 
to the cosmic rate production rate in the whole Universe, assuming the energy generation 
rate evolves rapidly with redshift following the luminosity density evolution of quasi-stellar 
sources. They obtain an upper limit which is ~ 3 Eqr. It is interesting however to mention 
that Ecr is comparable to that produced in 7-rays by cosmological GRBs (which are not 
expected to be absorbed by the intervening backgrounds since the typical photon energy is 
below 1 MeV). Assuming the efficiency with which electrons (which ultimately produce the 
observed photons by synchrotron radiation) and protons is the same inside the GRB fireball, 
the value of Eqr quoted above might well be closer to the actual value. For this, and other 
reasons, it is anyway unlikely that our calculations are accurate to better than a factor 3 or 
so. Moreover our neutrino event rate calculation might be conservative in this respect. 
The neutrino flux is given by 



dE v \ ^ ; E u > E B 

where A is a normalization constant that is determined from energy considerations as ex- 
plained above and in reference [0] . The observed neutrino rates are calculated by folding this 
generic flux over the distributions of individual bursts in distance, f(z), energy, f'(E GRB ) 
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and boost factor, f"(T): 

N„<xj J J J ^(E GRB , z, T, E v ) P(E U ) f(z) f'(E GRB ) f"(T) dE GRB dz dT dE v . (2) 

where E GRB is the energy emitted by a particular GRB, z its redshift and V the boost 
factor. P(E U ) represents the efficiency of detecting a neutrino of energy E v . The first two 
distributions can be modelled after observations: a cosmological distribution in distance, 
and an energy distribution which assumes that ten percent of GRB produce more energy 
than average by a factor of ten, and one percent by a factor of 100 0. 

Most important however are the fluctuations in the T factor around its average phe- 
nomeno logical value of 10 2 -10 3 . The fluctuations in T affect the efficiency for producing 
pions in the p — 7 collisions in the fireball as T -4 [II], as well as the break energy, Eb, which 
varies as T 2 . Unfortunately the distribution in T cannot even be guessed at. Nevertheless, 
it is critical in making quantitative predictions . The physics is clear. In GRBs, high 
luminosities are emitted over short times, therefore the large photon density would render 
GRB opaque unless T is very large. Only transparent sources with large boost factors emit 
photons. They are however relatively weak neutrino sources because the actual photon tar- 
get density in the fireball is diluted by the large Lorentz factor. An even moderately reduced 
value of T will produce a prolific neutrino source. 

We remind the reader that the results obtained from Eq. ^| are at variance with the neu- 
trino rate obtained by multiplying the average rate per burst by the number of bursts. Even 
neglecting the all important fluctuations in T, there is no such concept as an average GRB. 
E.g. for T = 300, the correct computation of Eq. |2| yields a rate of ~ 75 events per km 2 and 
year, roughly an order of magnitude larger than the prediction obtained by neglecting the 
observed burst-to-burst fluctuations in distance and energy. Another consequence of fluctu- 
ations is that the signal is dominated by a few very bright bursts, which greatly simplifies 
their detection. 

As pointed out above, the average neutrino energy varies with the square of boost factor 
and therefore the calculated event rates, especially their dependence on T, is strongly affected 
by the fact that higher energy neutrinos are preferentially absorbed in the Earth before 



reaching the detector Also this effect has been neglected in all previous calculations. For 
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instance, the 75 events just mentioned are reduced by a factor 3 by absorption. One should 
on the other hand remember that by oscillations, a large fraction of can oscillate into v T 



which penetrate the Earth |L^] . It is in this context important to realize that a kilometer-scale 
detector such as IceCube can measure the energy of the neutrinos. Therefore, signal events 
can be separated from the low energy atmospheric background by energy measurement, which 
allows the identifications of neutrinos from all directions and not just in the hemisphere where 
they pass through the Earth. 



3 Results 

We calculate the flux of neutrinos from GRB in the fireball model following reference [[F]. 
The number of protons in the fireball is fixed by the assumption that they are the source 
of the ultra high energy cosmic rays above 10 19 eV. The results are shown in Figure [1 
which shows the + flux from GRBs for different values of T. The fluxes have been 
multiplied by E 2 so that they represent the energy emitted in the form of neutrinos. Notice 
the variation of the break in the spectrum Eb as V 2 . Also, for values of V below ~ 100 the 
fireball becomes opaque to protons and at this point the total amount of energy available for 
neutrino production is converted. The neutrino flux, which roughly scales as T~ 4 , saturates 
and no longer grows with decreasing boost factor; see Figure |1| 

As usual we calculate the event rate per GRB by folding the neutrino flux with the 



probability of detecting a muon produced in a or interaction. Fluctuations in distance 
z and total GRB energy, or fluency, -Egrb are accounted for by a Monte Carlo simulation we 
have developed for this purpose. We simulate a large number of GRBs at different zenith 
angles assuming an isotropic distribution and subsequently obtain the event rate per year by 
assuming 1000 GRBs/year. Absorption of the neutrinos in the Earth prior to reaching the 
detector is taken into account using the density profile of reference ||13|| . It is also important 
to implement the fact that above the detector there is a limited column density of atmosphere 
and ice available for neutrino detection. 

In Table [l] we separately show, for different T's, the event rate of upgoing, i.e. neutrinos 
that cross the Earth before interacting near the detector, and downgoing neutrinos. The 
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E„ (GeV) Boost factor T 

Figure 1: Left: + z/ M fluxes from GRB for different values of the Lorentz factor T assuming 
GRBs are responsible for the observed cosmic ray spectrum above 10 19 eV. Right: Energy 
injected in GRB neutrinos as a function of the boost factor T of the fireball. 

first two columns show for comparison the number of events when neither absorption nor 
the limited amount of target above the detector are taken into account. It is clear that both 
effects play an important role in obtaining the correct event rate. This is not surprising: the 
Earth becomes opaque to neutrinos of energy around 100 TeV, and the muon range exceeds 
the ~ 2 km vertical depth of the IceCube detector at energies around 1 TeV. 

The dependence of the number of events on the T factor is shown in Figure 0. Two 
competing effects determine the shape of the curve. The event rate decreases with increasing 
T following the dependence of the neutrino flux which varies as T -4 . This decrease is partially 
offset because higher energy neutrinos resulting from larger boost factors are more efficiently 
detected. For low values of T, below about 100, the saturation of the total energy available 
for neutrino production is seen. On the other side of the T range, for large values of T, the 
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Events/ (km 2 yr) in 2n sr 


No absorption 


Absorption 


Downgoing 


Upgoing 


Downgoing 


Upgoing 


r = 100 


1133 


1112 


476 


600 


r = 300 


38 


38 


13 


14 


r = 1000 


0.14 


0.15 


4.2 x ID" 2 


2.8 x 10~ 2 



Table 1: + events (km 2 yr 1 ). Only fluctuations in distance and energy are taken into 
account. 

spectrum is very flat (~ E^ 1 ) up to ~ 70 PeV where the absorption by the Earth dominates. 
This reduces the event rate of upgoing events for values of V above 1000 as can be seen in 
Figure |2|. In the end the downgoing and upgoing event rates are similar except for large 
values of T. It is important to keep in mind that although downgoing neutrinos are not 
affected by absorption, their detection is limited by the column density of matter available 
for neutrino interaction. 

In Figure |3], we show the energy dependence of upgoing and downgoing + events 
for three representative values of T. The zenith angle distribution of upgoing neutrinos is 
shown in Figure |4], i.e. —1 < cos((9 zcn i t h) < 0. As T increases, the higher energy neutrinos are 
attenuated by the Earth. This explains why, as T increases, the distributions increasingly 
resemble an exponential attenuation function. 



3.1 v T + v T events 

v T production is expected to be very small in GRBs and in general in any astrophysical 
environment where z/s are produced in p — p or p — 7 collisions. Several calculations suggest 



a ratio 14 



Oscillation scenarios in which u^s convert into z/ r 's can however provide abundant sources 
of z/ r 's. Assuming a typical value F Vt / v = 0.5 suggested by SuperKamiokande measurements 
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Figure 2: + event rate per km 2 yr as a function of the boost factor T taking into 
account fluctuations in distance and the GRB fluency. Absorption in the Earth and the 
limited target above the detector are taken into account. 



n| , we obtain the double bang [JHJ [] event rates shown in Table |2|. The probability of 
detecting a v T induced double bang event in IceCube is typically two orders of magnitude 
smaller than the probability of observing a at energies around 10 PeV [14, F7|. This 
explains the small values of the event rates in Table 0. In the event rate calculation we have 
accounted for the energy loss of the i/ T 's when they propagate along the Earth's interior. This 



produces a pileup of events around 100 TeV, as pointed out in fll2| , reducing the number 
of upgoing double bang events with respect to the downgoing ones. This is due to the 
probability of ^-induced double bang detection which is limited to a broad peak between 
~ 1 PeV and ~ 100 PeV (outside of which it is negligible). This also explains the fact that 
1 Events in which two separated showers can be identified, one initiated by the struck nucleon and the 
other by the decay of the r produced in the v T charged current interaction. 
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Figure 3: Energy distribution of + events for T =100, 300 and 1000. Solid (dashed) 
lines correspond to upgoing (downgoing) neutrinos. Only fluctuations in distance and energy 
are accounted for. Absorption in the Earth and the limited target above the detector are 
taken into account. 

the energy distribution of the event rate peaks in the vicinity of 10 PeV (see Figure |5|). 

We also calculated the v T + v T events that would be detected by the appearance of a r 
which decays to // just below the detector. The event rates are also shown in Table pj. In this 
case the energy distribution of the events peaks around the energy at which the events pileup 
(100 TeV); see Figure [5]. The probability of detecting the \i is ~ 17 % of the probability of 
detecting it in a interaction due to the branching ratio of the r to \i decay channel. This 
accounts for a factor ~ 6 difference between v T — > r — > [i in Table ^] and the + rate in 
Table § 

The expected event rates generated by both mechanisms (double bang and r — ► \i decay) 
have very different and characteristic zenith angle distributions. These are shown in Figure 
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Figure 4: cos(6> zcn ith) distribution of v^+v^ events for T =100, 300 and 1000. cos(# zen i t h) = — 1 
corresponds to upgoing neutrinos and cos(6 , ze nith) = to horizontal. Notice that for T = 1000 
we have multiplied the event rate by a factor 1000. 

|5|. For large column depths inside the Earth (i.e. cos(6 l zenith ) ~ —1), the events pileup 
around 100 TeV and the small probability of detecting a double bang event at this energy 
reduces the number of double bangs with respect to events in which a /i is detected. When 
cos(# zen ith) ~ 0, because the amount of matter neutrinos have to cross is very small, there 
is no pileup of events and the angular distributions have roughly the same shape. Figure 
|6] also shows the zenith angle distribution of events produced by + for comparison, 
making it clear that despite of the flatter distribution of the /x's from v T + v T they are still 
outnumbered by /i's from interactions. 
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Events/ (km 2 yr) in 2tt sr 


Double Bang 


v T — > r — > /i 


Downgoing 


Upgoing 


Downgoing 


Upgoing 


r = 100 


0.54 


0.13 


38 


49 


r = 300 


3.6 x 10"2 


8.8 x 10~ 3 


1.2 


1.3 


r = 1000 


2.9 x 10" 4 


6.6 x 10~ 5 


3.6 x 10~ 3 


5.5 x lO" 3 



Table 2: v T + v T double bang events and events in which the r decays to /i (km~ 2 yr -1 ). 
Only fluctuations in distance and energy are considered. Absorption in the Earth and the 
limited column of matter above the detector are taken into account. 

3.2 Fluctuations in T 

Finally, we discuss burst to burst fluctuations in the boost factor T. So far we, quite un- 
realistically, assumed that T is single valued. We have, at present, no information on the 
distribution of T factors. It has been shown [0] that the neutrino rates can be significantly 
enhanced by fluctuations around the average value. It may be, of course, that the fluctua- 
tions in energy and distance which we took into account following the experimental evidence, 
already reflect some or all of the fluctuations in the boost factor. Following reference |7] we 
will illustrate the effect by assuming Gaussian distributions with half width o. The results 
are shown in Table |] for different values of (r). It is interesting to note that for (r) = 100 
the event rates are almost independent of the value of a. For T around 100 the event rate is 
weakly dependent on T due to the saturation of the amount of energy that goes into neutrino 
production (see Figures [l] and an d hence fluctuations in the value of the boost factor do 
not affect the event rate. This is not the case when (r) = 300 or 1000 since the event rate 
behaves roughly as T~ 4 in that T range. Figure |7] shows more clearly the dependence of the 
event rate on a for different values of (r) illustrating these points. It can be easily shown 
that when a <C (r) the following relation holds: 

Rate((r)±q) _a_ 
Rate((r)) K <r> ' U 

explaining why for (T) = 1000 the variation of the event rate with a is not as strong as when 

(T) = 300, even though the rate scales with T~ 4 in both cases. The same comments apply to 
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Figure 5: Energy distribution of upgoing i/ T + i? T double bang events (dashed line) and events 
in which the r decays to \l (solid line) for T =100, 300 and 1000. The double bang event 
rate has been multiplied by a factor 100 in all the plots. Absorption in the Earth and the 
limited target above the detector are taken into account. 

double bang events produced by i/ r 's and events in which the produced r decays to \i (both 
are shown in Table 

4 Summary and conclusions 

We have investigated potential signatures in km 3 telescopes of high energy neutrino fluxes 
produced in p — 7 interactions in GRB environments. We stress the fact that the rate is 
dominated by fluctuations in distance, GRB energy and in the bulk Lorentz factor T of the 
expanding GRB fireball. We have used an euclidean distribution to account for fluctuations 
in distance. Using a distribution in which GRB's follow the star formation rate, as suggested 
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Events/(km 2 yr) in 2n sr 


Absorption 


IDownffoinff 


I In Eroinff 


(r) = 100 


rr — n 
U — u 


4(0 


ouu 


a — 30 


47fi 


uw w 


a — 50 

yJ ',J\) 


472 


594 


a — 75 


464 


582 


rr — 1 00 

Ls X WW 


454 


571 


(r) = 300 


IT — 

\J — w 


lo 


1 /i 


a = 30 


15 


16 


a = 50 


22 


23 


a = 75 


39 


44 


cr = 100 


63 


75 


(r) = 1000 


a = 


4.2 x 10- 2 


2.8 x lO" 2 


a = 30 


4.2 x 10~ 2 


2.9 x 10~ 2 


a = 50 


4.3 x 10~ 2 


3 x lO" 2 


a = 75 


4.5 x 10~ 2 


3.1 x 10~ 2 


a = 100 


4.8 x 10~ 2 


3.4 x 10~ 2 



Table 3: + events (km~ 2 yr -1 ), taking into account fluctuations in boost factor T, 
distance and energy. The distribution of boost factors is assumed to be a Gaussian of half 
width a centered in (r) . Absorption in the Earth and the limited column of matter available 
for neutrino interaction above the detector are taken into account. 
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Events/(km 2 yr) in 2n sr 


Double Bang 


u T — > r — > fi 


Downgoing 


UDSoing 


Downgoing 


Upsoing 


(r) = 100 


(7 = 


0.54 


0.13 


38 


49 


a = 30 


0.52 


0.12 


37 


48 


a = 50 


0.48 


0.11 


36 


46 


a= 70 


0.43 


0.1 


33 


43 


a = 100 


0.39 


0.09 


32 


40 


(r) = 300 


a = 


3.6 x 10~ 2 


8.8 x 10~ 3 


1.2 


1.3 


a = 30 


4.0 x 10~ 2 


1.0 x 10~ 2 


1.4 


1.5 


a = 50 


5.0 x 10~ 2 


1.2 x 10~ 2 


1.9 


2.1 


a = 70 


6.9 x 10~ 2 


1.6 x 10~ 2 


2.9 


3.3 


a = 100 


9.9 x 10~ 2 


2.4 x 10~ 2 


5.0 


6.2 



Table 4: v T + v T double bang events and events in which the r decays producing a ji 
(km -2 yr -1 ). We take into account fluctuations in boost factor T, distance and energy. The 
distribution of boost factors is assumed to be a Gaussian of half width a centered in (r). 
Absorption in the Earth and the limited column of matter available for neutrino interaction 
above the detector are taken into account. 
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Figure 6: cos(^ zen ith) distribution of u T + v T double bang events (dashed line) as well as 
events in which the r decays to \i (solid line), for T =100, 300 and 1000. cos(# zcn ith) = 
— 1 corresponds to upgoing neutrinos and cos(6 , zonit h) = to horizontal. Also shown for 
comparison are the + event rates (dotted line). 



in | 20fl , the event rates are reduced only by ~ 20%. On the other hand, a cosmological 



distribution following galaxies |L8], [T9|] allows for more close GRB's so that the event rate 
increases by a factor of 20 (see Table 0). Neutrino telescopes may help to constrain the 
distribution in distance of GRBs. The most relevant parameter is T, which determines the 
rate of p — 7 interactions and hence the amount of energy that goes into neutrino production. 
For small values of T (T ~ 50 or less) the expansion of the GRB fireball is not sufficiently 
fast and the large photon density makes it opaque to p — 7, efficiently producing pions - 
the parents of the neutrinos — and saturating the amount of energy available for neutrino 
production. Due to the scaling of the energy break in the spectrum with T 2 , mostly low 
energy neutrinos are produced whose detection efficiency is smaller due to the small range of 
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Figure 7: Upgoing v^+Vp event rate as a function of a, the half width of the boost distribution 
around (r), for (r) = 100, 300 and 1000. The rate for (T) = 1000 has been scaled up by a 
factor 100. Absorption in the Earth is taken into account. 

the muon at low energies and the relatively high energy threshold of the neutrino telescopes 
(~ 100 GeV). However the large neutrino flux compensates for the small detection efficiency 
For large values of T (T > 300), an increasingly larger fraction of the neutrino energy goes 
into the high energy part of the spectrum, however the overall amount of energy is very 
small, producing small event rates. Values of T < 100 give rise to large number of events 



that would even be observable in existing smaller neutrino telescopes such as AMANDA |2T| 
(one has to scale the results down by roughly two orders of magnitude to account for the 
smaller effective area of the detector). 

We have shown that absorption of the upgoing u^s inside the Earth, as well as the limited 
column of matter available for downgoing neutrino interactions, play a relevant role, making 
upgoing and downgoing event rates roughly equal. The change of the energy break of the 
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Events/(km 2 yr) in 2n sr 


Euclidean 


Cosmological 


Star Formation Rate 


r = 100 


1076 


21,029 


832 


r = 300 


27 


424 


20 


r = 1000 


0.3 


1.8 


6 x 10~ 2 



Table 5: -\- events (km -2 yr -1 ). Only fluctuations in distance and energy are taken 
into account. Different distributions in distance are used: the first column corresponds to 
an euclidean distribution, the second to a cosmological distribution following galaxies ]18[ 



and the third one assumes that GRB's follow the distribution of star formation regions [^0 
Absorption in the Earth as well as the limited column of matter above the detector are taken 
into account. 

spectrum with T 2 combined with the absorption of the Earth is reflected in the zenith angle 
distributions of the event rates which may give some complementary information about T. 

GRB neutrino detection with km 3 neutrino telescopes also has the potential to investigate 
v t oscillations over a baseline of 1000 Mpc. Double bang v T induced events offer an 
unmistakeable signature which allows downgoing v T detection. A km 3 telescope operating 
for 10 years may detect ~ 10 downgoing double bang events if T = 100 without any potential 
background. Upgoing double bang events are not going to be detected since the iv's pileup 
around 100 TeV, where the probability of a double bang detection is negligible. Muons 
produced in r decays are outnumbered by /i's from interactions at least for the type of 
fluxes expected from GRBs |p7| ; besides, their energy distribution does not show a clear and 
characteristic signature so they are difficult to identify. 

In summary, neutrino telescopes open up the possibility of determining the value of T and 
its fluctuations, as well as the possibility of identifying iv's. They are powerful instruments 
to reveal important astrophysical information about the most energetic objects ever observed 
in the universe and about neutrino oscillation scenarios over cosmological baselines. 
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